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Abstract 

^ ' A novel mechanism for sympatric speciation that takes into ac- 

j^ ■ count complex bio-processes within each individual organism is pro- 

posed. According to dynamical systems theory, organisms with iden- 
tical genotypes can possess differentiated physiological states and may 
coexist 'symbiotically' through appropriate mutual interaction. With 
mutations, the phenotypically differentiated organisms gradually come 
to possess distinct genotypes, while maintaining their symbiotic rela- 
tionship. This symbiotic speciation is robust against sexual recombi- 
nation, because offspring of mixed parentage, with intermediate geno- 
types, are less fit than their parents. This leads to sterility of the 
hybrid. Accordingly, a basis for mating preference also arises. 



The question posed by Darwin(1859) of why organisms are separated into 
distinct groups, rather than exhibiting a continuous range of characteristics, 
has not yet been fully answered. In spite of several explanations involving 
sympatric speciation, "we are not aware of any explicit model demonstrat- 
ing the instability of a sexual continuum," according to Maynard-Smith and 
Szathmary(1995). The difficulty involving stable sympatric speciation is that 
it is not clear how two groups, which have just begun to separate, coexist 
while mutually interacting. Here, we propose a mechanism through which 
two groups with little (or no) difference in genotype form a 'symbiotic' rela- 
tionship under competition. This mechanism is understood in terms of the 
'isologous diversiffcation theory' (Kaneko & Yomo 1997,1999), according to 
which organisms with identical genotypes spontaneously split into distinct 
phenotypes and establish a relationship in which the existence of one group 
is mutually supported by the other. By considering genetic mutations and 
sexual recombinations, a sympatric speciation process follows, resulting in 
the formation of distinct genotypic groups exhibiting reproductive isolation. 
This process is robust with respect to ffuctuations, due to the symbiotic 
relationship. The hybrid offspring of the two groups becomes sterile, and 
also provides a basis for mating preference, a major mechanism in sympatric 
speciation(Maynard-Smith 1966, Lande 1981, Turner & Burrows 1995, Kon- 
drashov & Kondrashov, 1999, Dieckmann & Doebeli 1999, Futsuyma 1986, 
Howard & Berlocher ed. 1998). 

To study phenotypic and genotypic diversification through interaction, 
we have to consider a developmental process that maps a genotype to a phe- 
notype. As an illustrative model, we consider a dynamic process consisting of 
several interacting metabolic cycles. Each organism possesses such internal 
dynamics with several metabolic cycles, while it selectively consumes exter- 
nal resources, depending on their internal dynamics, and transforms them 
into some products. Through this process, organisms mature and eventually 
become ready for reproduction. 

For most studies on population biology and evolution so far, it is widely 
assumed that a phenotype is uniquely determined for a given genotype and 
environment. If this assumption were always true, population dynamics only 
of genotypes would be sufficient to study the evolutionary process theoret- 
ically. However, there are cases that organisms of the same genotype may 
take distinct phenotypes through interaction. 

First, some mutant genotypes related to malfunctions show various phe- 



notypes, each of which appears at a low probabihty(Holmes 1979). This 
phenomenon is known as low or incomplete penetrance (Opitz 1981), which 
suggests plastic ontogenesis. 

Although the origin of low penetrance in multicellular organisms may 
not be well clarified, differentiation of physiological states is already known 
in bacteria (Novick and Weiner, 1957). Furthermore, one of the authors and 
his colleagues have found that specific mutants of E. coli show (at least) two 
distinct types of enzyme activity, although they have identical genes. These 
different types coexist in an unstructured environment of a chemostat (Ko et 
al. 1994), and this coexistence is not due to spatial localization. Coexistence 
of each type is supported by each other. Indeed, when one type of E. coli is 
removed externally, the remained type starts differentiation again to recover 
the coexistence of the two types. In addition, even at a molecular level, 
a mutant gene of xylanase was shown to produce various levels of enzyme 
activity(Ko et al., 1994). A mechanism for a single gene to show various 
levels of molecular function has also been elucidated in physico chemical terms 
(Kobayashi et al. 1997). 

Such differentiation of phenotype is also discussed as a possibility of dif- 
ferent inheritable states of a same genotype (see e.g.. Landman 1991). Al- 
though we do not assume any epigenetic inheritance here since its relevance 
to evolution is still controversial, it should be noticed that the existence of 
different physiological states from a single genotype itself is demonstrated 
experimentally, even if one does not accept the inheritance of the states. 

Note also that the differentiation of phenotypes from the same genotype 
is supported also theoretically, as will be mentioned later. Here we will study 
the relevance of such phenotypic diversification to evolution. 

Model 

In our theoretical model, the phenotype is represented by a set of vari- 
ables, corresponding to metabolic processes or some other biological pro- 
cesses. To be specific, each individual i has several (metabolic or other) cyclic 
processes, and the state of the j-th process at time n is given by Xl{i). With 
k such processes, the state of an individual is given by the set (X^(i),X^(i), 
■ ■ ■ ,X!^{i)), which defines the phenotype. This set of variables can be re- 
garded as concentrations of chemicals, rates of metabolic processes, or some 
quantity corresponding to a higher function. The state changes temporally 
according to a set of deterministic equations with some parameters. 

Since genes are simply chemicals contained in DNA, they could in prin- 



ciple be included in the set of variables. However, according to the 'central 
dogma of molecular biology' (Alberts et al. 1994), the gene has a special role 
among such variables: Genes can affect phenotypes, the set of variables, but 
the phenotypes cannot directly change the code of genes. During one gener- 
ation, changes in genes are negligible compared with those of the phenotypic 
variables they control. Hence, in our model, the set corresponding to genes is 
represented by parameters that govern the dynamics of phenotypes, since the 
parameters in an equation are not changed while they control the dynamics 
of the variables. 

Our model consists of the following procedures. 

(i) Dynamics of the phenotypic state: The dynamics of the variables Xl{i) 
consist of a mutual influence of cyclic processes (X^(i)) and interaction with 
other organisms (X^(i')). 

(ii) Growth and Death: Each individual splits into two when a given 
condition for growth is satisfied. Taking into account that the cyclic pro- 
cess corresponds to a metabolic, genetic or other process that is required 
for reproduction, we assume that the unit replicates when the accumulated 
number of cyclic processes goes beyond some threshold. To introduce com- 
petition for survival, death is included both by random removal of organisms 
at some rate as well as by a given death condition based on their state. 

(iii) Genetic parameter and mutation: Following the above argument, 
genes are represented as parameters in the model. With reproduction, the 
genes slowly mutate. The set of parameters in the model changes slightly 
through mutation when offspring is reproduced. 

To be specific we consider the following model. First, the state variable 
Xl{i) is split into its integer part -R^(i) and the fractional part a:^(i) = 
mod[X!^{i)]. The integer part Ri{i) is assumed to give the number of times 
the cyclic process has occurred since the individual's birth, while the frac- 
tional part x^^{i) gives the phase of oscillation in the process. As a simple 
example, the internal dynamics of the cyclic process is assumed to be given 
by J2m^^sin(27ix^{i)), while the interaction among organisms is given by 
the competition for resources among the Nn organisms existing at the mo- 

ment, given by 1 \i) = psin[27ix^[i)) H ^—^ . (The second term 

comes from the constraint J2i -^^(^) = s^, i-e., the condition that A^ individuals 
compete for a given resource s^ at each time step. The first term represents 
the ability to secure the resource, depending on the state.) Our model is 



given by 

K+A<) = K(i) + E.„ ^sm(2,r<(,:)) - E™ !^sm{2rxi(i)) 

Then, as a specific example, the condition for the reproduction is given 
by '}liRn{i) > Thr. The rotation number -R^(i) is reset to zero when the 
corresponding individual splits. On the other hand, as the death condition, 
an individual with -R^(i) < — 10 (i.e., with a reverse process) is removed. 

Next, genotypes are given by a set of parameters a'"^(i), representing 
the relationship between the two cyclic processes d, and m {1 < i,m < k). 
With each division, the parameters a™^ are changed to a™^ + 6 with 6, a 
random number over [— e, e], with small e, corresponding to the mutation 
rate. Although the results in the figure adopt the mutation rate e = .001, 
change of mutation rate is just responsible for the speed of the separation of 
the parameters, and the conclusions are independent of its specific value. 

Let us make some remarks on our model. Each term a^™sin(a;™(i)) gives 
how a process x™ influences on x^. For example, in a metabolic process, one 
cycle influences some other through catalytic reactions, depending on the 
activity of an enzyme corresponding to it. With the change of genes, the 
activities of enzymes can change, which leads to the change of the parameter 
values of a^™ accordingly. Following this argument, genotypes are regarded 
to be represented by a set of parameters a^"^. Indeed, we have also studied a 
specific biochemical reaction network model with its catalytic efficiency as a 
genetic parameter, and the results to be presented are observed. 

The interaction term psin{2Tixl^{i)) represents the influence on the cyclic 
process between individuals through the exchange of chemicals (or by other 
means). Since this term can change its sign, chemicals can be secreted from 
each individual to the environment. Then, the resources that are taken out 
from one individual may be used by some other. Through this ecological 
interaction, individuals may keep some relationship if they are differentiated. 

Of course, the above explanation is just one example of correspondence 
of our model to a real biological process. As long as its mathematic structure 
is common, the validity of our model is not restricted to the above example, 
and the results to be presented can generally be applied. 

Scenario for Symbiotic Sympatric Speciation 



The above model is one of the simplest to discuss loose developmental 
process. We have also carried out simulations of several models of this type, 
for example, consisting of a metabolic process of autocatalytic networks. 
Through the simulations and theoretical considerations of several models 
the following mechanism for speciation is proposed. Since the characteristic 
features for speciation, to be presented as follows, are common, we adopt 
the above simplest model to illustrate the scenario here. Note, of course, 
the scenario for the speciation is expected to work in a more realistic model 
including much complicated processes for development and interaction. 

Stage-1: Interaction-induced phenotypic differentiation 

When there are many individuals interacting for finite resources, the phe- 
notypic dynamics begin to differentiate even though the genotypes are iden- 
tical or differ only slightly (see the 'light blue points' in Fig.l. This differenti- 
ation generally appears if nonlinearity is involved in the internal dynamics of 
some phenotypic variables(Kaneko & Yomo 1997,1999, Furusawa & Kaneko 
1998). Then, slight phenotypic differences between individuals are amplified 
by the internal dynamics (e.g., metabolic reaction dynamics). Through in- 
teraction between organisms, the different phenotypic dynamics tend to be 
grouped into two (or more) types, despite the fact that all have identical (or 
only slightly different) genotypes. In the example of Fig. la, the phenotype 
splits into two groups, which we refer to as 'upper' and 'lower' groups. 

This differentiation process has recently been clarified as isologous di- 
versification(Kaneko & Yomo 1997,1999), in which two groups with distinct 
phenotypes appear even in a group with a single genotype. This interaction- 
induced phenotype diversification is a general consequence when the develop- 
mental process with interaction between organisms is considered as a nonlin- 
ear dynamics process(Kaneko 1991,1994). When there is an instability in the 
dynamics, the temporal evolution of individuals in phenotype space begins to 
diverge. Then, through interactions, these dynamics are stabilized through 
the formation of distinct groups with differentiated states in the pheno-space. 
The existence of the two (or more) groups eliminates the instability in the 
dynamic (metabolic) process that exists when one of the groups is isolated. 
Hence, the existence of all groups is required for the survival of each. For 
example, if a group of one type is removed, then the phenotype of individuals 
of another type changes in compensation for the missing type. 

To put the above explanation in biological terms, consider a given group 
of organisms faced with a new environment and not yet specialized for the 
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processing of certain specific resources. Each organism has metabohc (or 
other ) processes with a biochemical network. As the number of organisms 
increases, they compete for resources. The interaction, for example, results 
from the use of some byproduct of one organism by others. As this interaction 
becomes stronger, the phenotypes become diversified to allow for different 
uses of metabolic cycles, and they split into two (or several) groups. Each 
group is specialized in some metabolic cycles, and also in the processing of 
some resources. Here, the byproduct of metabolic processes of one group is 
necessary for another group to allow for it to specialize in some particular 
metabolic cycle. Resources secreted out from one group can be used as a 
resource for the other group, and vice versa. Hence, the two groups realize 
a differentiation of roles and form a symbiotic relationship. Each group is 
regarded as specialized in a different niche, which is provided by another 
group. 

As an extreme case, this differentiation can occur even when a single 
resource is supplied externally (i.e., s^ = for 2 < j < A; = 3). In this case, 
the temporal average oi psin{x'^{i)) is positive for one group and negative 
for the other, while that of psin{x\{i)) has the opposite sign. With this 
differentiation of phenotypes, one group uses x^ as a resource for growth 
provided by the other, which uses x^ as a resource. 

It should be pointed out that the the progeny of a reproducing individual 
belonging to one group may belong to the other group at this stage, since all 
groups still have almost identical gene sets. 

Stage-2: Co-evolution of the two groups to amplify the differ- 
ence of genotypes 

Now we discuss the evolutionary process of genotypes. After the phe- 
notype is differentiated into two groups, the genotype (parameter) of each 
group begins to evolve in a different direction, as shown in Fig.l and Fig.2. 
This evolution occurs, since for the upper (lower) group, those individuals 
with a smaller (larger) parameter value reproduce faster. In our numerical 
simulations, there always exists (at least one) such parameter. As a simple 
illustration, assume that the two groups use certain metabolic processes dif- 
ferently. If the upper group uses one metabolic cycle more, then a mutational 
change of the relevant parameter to enhance this cycle is favored for the up- 
per group, while a change to reduce it (and enhance some others) may be 
favored for the lower group. In other words, each organism begins to adapt 
in one of the niches formed by another (or others). 
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Note that R^ also takes a different value between the two groups in the 
other way round, since J2j R^ = Thr for each division. As for the parameter 
change, the values of a^^ and a^^ split first in this example, but then a^^ and 
later other parameter values also start to split into the two groups. Several 
genes start to be responsible for the differentiation. 

As the genetic separation progresses, phenotypic differences between the 
two groups also become amplified (see Fig.l). With the increase in the split 
in genotypes, it begins to become the case that offspring of members of a 
given group certainly keep the phenotype of this group. Since the phenotype 
of one group stabilizes the other, the evolutions of the two groups are inter- 
dependent. Hence the tempo of the genetic evolution in one group is related 
to that of the other. The two groups co-evolve, maintaining their 'symbiotic' 
relationship, established in the previous stage. Indeed, as shown in Fig. 2, 
the growth speeds of the two groups remain of the same order, even if each 
genotype and phenotype change with time. 

With the above described co-evolutionary process, the phenotype differ- 
entiation is fixed to the genotype. In much later generations, this fixation is 
complete. In this case, even if one group is isolated, offspring with the phe- 
notype of the other group are no longer produced. Offspring of each group 
keep their phenotype (and genotype) on their own. 

Reproductive Isolation with respect to Sexual Recombination 

The importance of the present scenario lies in the robustness of the spe- 
ciation process. Even if one group happens to disappear through some fluc- 
tuations at the initial stage of the speciation process, coexistence of the two 
distinct phenotypic groups is recovered. Hence, the present process is also 
expected to be stable against sexual recombination, which mixes the two 
genotypes and may bring about a hybrid between the two genotypes. To 
demonstrate this stability, we have extended the previous model to include 
this mixing of genotypes by sexual recombination. 

Here, we have modified our model so that the sexual recombination occurs 
to mix genes. To be specific, the reproduction occurs when two individuals ii 
and ^2 satisfy the threshold condition (J^e^niik) > Thr), and then the two 
genotypes are mixed. As an example we have produced two offspring j = ji 
and J2, from the individuals ii and ^2 as 

a^™(j) = a^™(ii)r^^™ + a^"(i2)(l - rj™) + 6 (2) 



with a random number < r!^™ < 1 to mix the parents' genotypes, besides 







the random mutation term by 5. Even if two separated groups may start to 
be formed according to our scenario, the above recombination forms 'hybrid' 
offspring with intermediate parameter values a^™ when two organisms from 
different groups mate. Also, depending on the random number, for some 
offspring, the parameter value a'™ may be closer to one of the parents, but 
that of a} ™ may be closer to that of the other of the parents. Accordingly 
recombinations of the two group can lead to a different combination of alle- 
les, since the two groups take different combination of the parameter values 
{a^™}, from the two groups. 

Although the hybrid is formed in this random mating with some propor- 
tion (1/2 if the two groups have equal population), it turns out that this 
hybrid, irrespective of which phenotype it realizes, has a lower reproduction 
rate than the other two groups which have a 'matched' genotype-phenotype 
correspondence with a higher reproduction rate. In Fig. 3, we have plotted 
the average offspring number for given genotype parameters (to be precise, 
the average over a given range of parameters). As shown, a drop at the inter- 
mediate value in the offspring number starts to appear, through generations. 
Within few dozen generations, as certain genotypic parameters are apart, 
there is little or no chance for a hybrid to reproduce, and Fl sterility results. 

Note that this conclusion is drawn even without assuming mating pref- 
erence. Rather, it is natural, according to the present scenario, that mating 
preference in favor of similar phenotypes evolves, since it is disadvantageous 
for individuals to produce a sterile hybrid. In other words, the present mecha- 
nism also provides a basis for the evolution of sexual isolation through mating 
preference. Note, however, that in sympatric speciation starting from only 
the mating preference, one of the groups may disappear due to fluctuations 
when its population is not sufficiently large. In contrast, according to our sce- 
nario, the coexistence of the two groups is restored even under disturbances. 
Hence, it is concluded that our mechanism yields robust sympatric specia- 
tion, i.e., differentiation of geno- and pheno- types and sexual reproductive 
isolation(Dobzhansky, 1951), even in the situation in which all individuals 
interact with all others equally. 

Importance of Phenotypic Differentiation 

Evolution according to our scenario often leads to specialization with re- 
gard to resources through competition. Indeed, the coexistence of two (or 
more) species after the completion of the speciation is also supported by 

9 



the resource competition theory of Tilman(1976,1981). However, to reahze 
the speciation process, phenotype differentiation from a single genotype is 
essential. As long as phenotype is uniquely determined by genotype, two in- 
dividuals with a slight genotype difference can have only a slight phenotype 
difference also. Since competition is strong among individuals with simi- 
lar phenotypes, they cannot coexist as a different group. Hence two groups 
cannot be differentiated from a group of single (or similar) genotypes. Con- 
trastingly, in our scenario, even if the genotypes of two individuals are the 
same or only slightly different, their phenotypes need not be similar, and can 
in fact be of quite different types, as shown in Fig.l. Accordingly, these two 
groups can coexist. 

To check the importance of this phenotypic differentiation from a single 
genotype, we have also performed several numerical experiments with our 
model, by choosing parameters so that differentiation into two distinct phe- 
notype groups does not occur initially. In this case, separation into two (or 
more) groups with distinct pheno/geno-types is not observed, even if the ini- 
tial variance of genotypes is large, or even if a large mutation rate is adopted. 
This clearly demonstrates the relevance of phenotypic differentiation. 

On the other hand, the genetic differentiation always occurs when the 
phenotype differentiates into two (or more) groups. To be specific, in our 
model, the condition for the differentiation is as follows: First, the param- 
eter p should be larger than some value. For example, for A; = 3 with 
s^ = 2,s^ = 4,s^ = 6 and with the initial parameters a^"^{i) ~ — 0.2/(27r), 

the differentiation appears for p ^ 1.8. Second, the resource term per unit 
{J2j s^ /N) should be smaller than some threshold value. For example, the 
threshold resource is Sthr ~ 10, for 3^ = 3"^ = 3^, p = 1.5/(27r), A^ ^ 300 and 
the initial parameters a^"^{i) ~ — .l/(27r). Note that these conditions imply 
strong interaction in competing for resources, and are easier to be satisfied, 
as the number of individuals competing for given resources increases. 

Discussion 

In the symbiotic speciation process, the potentiality for a single genotype 
to produce several phenotypes declines. After the phenotypic diversifica- 
tion of a single genotype, each genotype again appears through mutation 
and assumes one of the diversified phenotypes in the population. Thus the 
one-to-many correspondence between the original genotype and phenotypes 
eventually ceases to exist. As a result, one may see a single genotype ex- 
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pressing small numbers of phenotypes in nature, since most organisms at the 
present time have gone through several speciation processes. One can also 
expect that mutant genotypes tend to have a higher potentiality than the 
wild-type genotype to produce various phenotypes. Indeed, this expectation 
is consistent with the observation that low or incomplete penetrance(Opitz 
1981, Holmes 1979) is more frequently observed in mutants than in a wild 
type. 

Taking our results and experimental facts into account, one can predict 
that new species or organisms emerging as a species have a high potentiality 
to produce a variety of phenotypes, while 'living fossils', such as Latimeria 
chalumnae and Limulus, have a stable expression of a small number of phe- 
notypes. Relationship between evolvability and plasticity in ontogenesis is 
an important topic to be pursued. 

Since the speciation discussed in this paper is triggered by interaction, and 
not merely by mutation, the process is not so much random as deterministic. 
In fact, the speciation process occurs irrespectively of the adopted random 
number in the simulation. Some of the phenotypic explosions in nature that 
have been recorded as occurring within short geologic periods, may have 
followed the deterministic and relatively fast process of interaction-induced 
speciation. Hence, our scenario may shed some light on the variation of 
timescales on which evolution proceeds, e.g., punctuated equilibrium(Gould 
& Eldredge 1977). Here it should be noted that the change in phenotypes 
occur in few generations. The speed of genetic change, of course, depends on 
the mutation rate, but the present mechanism is found to work even for any 
smaller mutation rate (say e = 10~^). 

In the present paper, we have mostly reported the case only with 3 process 
{k = 3), but we have numerically confirmed that the present speciation 
process works also for k > 3 (e.g., k = 10). By choosing a model with 
many cyclic processes, we have also found successive speciation of genotypes 
into several groups from a single genotype. With evolution, the phenotypes 
begin to be separated into two groups, each of which is specialized in some 
processes, and depends on the byproducts of the other. Later, the species 
diverge into further specialized groups, which are fixed into genotypes. This 
process is relevant to consider adaptive radiation. 

Discussion of the mechanism involved in evolution often consists of mere 
speculation. Most important in our scenario, in contrast, is its experimen- 
tal verifiability. Isologous diversification has already been observed in the 
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differentiation of enzyme activity of E. coli with identical genes(Ko et al. 
1994). By observing the evolution of E. coli in the laboratory(Xu et al., 1996, 
Kashiwagi et al., 2000), controlling the strength of the interaction through 
the population density, one can check if the evolution on genetic level is 
accelerated through interaction-induced phenotypic diversification. Our isol- 
ogous symbiotic speciation, based on dynamical systems theory, numerically 
confirmed and biologically plausible, can be verified experimentally. 

Acknowledgment: We thank M. Shimada for illuminating comments. The 
present work is supported by Grants-in-Aids for Scientific Research from the 
Ministry of Education, Science and Culture of Japan (Komaba Complex 
Systems Life Project). 



References 

• B. Alberts et al.. The Molecular Biology of the Ceii(Garland, New York 
ed.3. 1994). 

• C.Darwin, On the Origin of Species by means of natural selection or the 
preservation of favored races in the struggle for life (Murray,London,1859). 

• U.Dieckmann and M.Doebeli, On the origin of species by sympatric 
speciation, Nature400, 354 (1999). 

• T.Dobzhansky, Genetics and the Origin of Species (Columbia Univ. 
Press., New York, ed.2, 1951). 

• C.Furusawa and K.Kaneko, Emergence of Rules in Cell Society: Dif- 
ferentiation, Hierarchy, and Stability, Bull. Math. Biol.QO, 659(1998). 

• D. J.Futsuyma, Evolutionary Biology (Sinauer Associates Inc., Sunder- 
land, Mass, ed.2, 1986). 

• S.J.Gould and N.Eldredge, Punctuated equilibria: the tempo and mode 
of evolution reconsidered. Paleobiology 3, 115 (1977). 

• L.B.Holmes, Penetrance and expressivity of limb malformations. Birth 
Defects. Orig. Artie. Ser. 15, 321(1979). 

12 



D.J. Howard and S.H.Berlocher,Ed., Endless Form: Species and Speci- 
ation (Oxford Univ. Press., 1998). 

K.Kaneko, Clustering, Coding, Switching, Hierarchical Ordering, and 
Control in Network of Chaotic Elements, Pi2ysica41D, 137(1990). 

K.Kaneko, Relevance of Clustering to Biological Networks, Piiysica75D,55(1994). 

K.Kaneko and T.Yomo, Isologous Diversification: A Theory of Cell 
Differentiation, BulLMath.Biol.59,139{1997). 

K.Kaneko and T.Yomo, Isologous Diversification for Robust Develop- 
ment of Cell Society, J. Theor. Bioi.l99,243(1999). 

A.Kashiwagi, T.Kanaya, T.Yomo, I. Urabe, How small can the dif- 
ference among competitors be for coexistence to occur. Researches on 
Population Ecology 40, in press. 

E.Ko, T.Yomo, I.Urabe, Dynamic Clustering of bacterial population, 
Physica 75D, 81 (1994) 

C. Kobayashi et al.. Thermal conversion from low- to high- activity 
forms of catalase I from Bacillus stearothermophilus J. Biol. Cheni. 
272, 23011 (1997). 

A.S.Kondrashov and A.F.Kondrashov, Interactions among quantitative 
traits in the course of sympatric speciation. Nature 400,351(1999). 

R.Lande, Models of speciation by sexual selection on phylogenic traits, 
Proc. Natl. Acad. Sci. USA 78, 3721 (1981). 

O.E. Landman, The inheritance of acquired characteristics, Annu. Rev. 
Genet. 25, 1 (1991) 

J.Maynard-Smith, Sympatric Speciation, The American Naturalist 100,637(1966). 

J.Maynard-Smith and E.Szathmary, The Major Transitions in Evolu- 
tion (W.H. Freeman, 1995). 

A. Novick and M. Weiner, Enzyme induction as an all-or-none phe- 
nomenon. Biochemistry 43, 553 (1957). 

13 



J.M.Opitz, Some comments on penetrance and related subjects, Am- 
J-Med-Genet. 8, 265 (1981). 

D.Tilman, Ecological competition between algae: Experimental confir- 
mation of resource-based competition theory, Science 192, 463 (1976). 

D.Tilman, Test of resource competition theory using four species of 
lake Michigan algae. Ecology 62, 802(1981). 

G.F.Turner and M.T. Burrows, A model for sympatric speciation by 
sexual selection, Proc. R. Soc. London B 260, 287 (1995). 

W.-Z.Xu, A.Kashiwagi, T.Yomo, I.Urabe, Fate of a mutant emerging 
at the initial stage of evolution. Researches on Population Ecology 38, 
231(1996). 



14 



o + n X <] 

oooo o 
oooo o 
oooo Ln 

1- -I- COLO (M 

^tV^io 
oooo 

OOOLO 
OOOCNi 
tDCDi- 
tDCO-i- 





o 
o 



o 
o 



IH 



Figure 1: Evolution of genotype-phenotype relationship. In the present 
model, due to the nonlinear nature of the dynamics, x^ often oscillates in 
time chaotically or periodically. Hence it is natural to use the integer part 
R^{j), as a representation of the phenotype, since it represents the number 
of cyclic process used for reproduction. Here (i?^,a^^) is plotted for every 
division of individuals. The first 2500 divisions are plotted in light blue, 
divisions 2501-5000 in pink, 11000-16000 in red, 36000-41000 in blue, and 
66000-71000 in green. Initially, phenotypes are separated, even though the 
genotypes are identical (or only slightly differ), as shown in light blue. Later, 
the genotypes are also separated, according to the difference in phenotypes. 
In the simulation, the population size fluctuates around 300, after an initial 
transient. (Hence the generation number is given roughly by dividing this 
division number by 300.) In the figures of the model, we adopt the following 
parameter values and initial conditions. The threshold number Thr for the 
reproduction is 1000, and the mutation rate of the parameters e is 0.001. 
Tnitiallv. the p-enotvne narameter.s are set as n^^ = —0.1 /f27rV The narameter 
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Figure 2: The evolution of the genotypic parameter. The parameter a^^(i) 
is plotted as a dot at every division (reproduction) event in (a), with the 
abscissa as the division number. The average time necessary for division 
(reproduction) is plotted for the upper and lower groups, where the average 
is taken over 2000 division events (6th - 8th generation). As the two groups 
are formed around the 2000th division event, the population size becomes 
twice the initial, and each division tinie is also approximately doubled. Note 
that the two average division speeds of the two groups remain of the same 
order, even when the genetic parameter evolves in time. 
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Figure 3: The average offspring number before death is plotted as a function 
of the parameter (genotype), for simulations with sexual recombination. As 
an extension to include sexual recombination, we have also studied a model 
in which two organisms satisfying the above threshold condition mate to 
reproduce two offspring. When they mate, the offspring have parameter 
values that are randomly weighted average of those of the parents, as given 
in the text. We have measured the number of offspring for each individual 
during its lifespan. By taking a bin width 0.005 for the genotype parameter 
a^^, the average offspring number over a given time span is measured to give a 
histogram. The histogram over the first 7500 divisions (about 20 generations) 
is plotted by the solid line (I), and the histogram for later divisions is overlaid 
with a different line, as given by II (over 7500-15000 divisions). III (1.5-2.25 
xlO^), IV(2.25-3 xlO^), and V(3.75-4.5 xlO^). As shown, a hybrid offspring 
will be sterile after some generations. Here we have used the same model 



